NMR Studies of pseudogap and electronic inhomogeneity in Bi2Sr2CaCu208+5 
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We report ^''O NMR measurements in single crystals of overdoped Bi2Sr2CaCu20g+i with Tc — 82 
K. We measure the full anisotropy of the planar oxygen Knight shift, electric field gradient, and spin 
lattice relaxation rate tensors, and show that the entire temperature dependence is determined by 
the suppression of the density of states in the pseudogap below T* --^ 94 K. The linewidth can be ex- 
plained by a combination of magnetic and quadrupolar broadening as a result of an inhomogeneous 
distribution of local hole concentrations that is consistent with scanning tunneling microscopy mea- 
surements. This distribution is temperature independent, does not break d symmetry, and exhibits 
no change below T* or Tc. 

PACS numbers: 74.25.nj, 74.72.-h, 74.81.-g 



I. INTRODUCTION 

The normal state of the high temperature supercon- 
ductors continues to attract attention despite more than 
twenty years of research. The parent state of these ma- 
terials is a Mott insulator, but for sufficient hole dop- 
ing the conductivity becomes metallic and high tem- 
perature superconductivity emerges with Tc ^ 100 K 
and d-wave symmetry. ^'■^ Other strongly correlated elec- 
tron systems such as the heavy fermions, the iron- 
based superconductors and the organic superconduc- 
tors also exhibit superconductivity in close proximity to 
an antiferromagnetism,''^^ but the cuprates are unique 
in that they also exhibit a pseudogap in the normal 
state over a broad range of dopings.^ It remains unclear 
whether the partial suppression of the low energy den- 
sity of states in the pseudogap is the result of a new 
thermodynamic phase^^^ or a crossover in magnetic be- 
havior driven by the proximate Mott insulating state. 
Of particular interest is the fact that the electronic de- 
grees of freedom in the pseudogap develop short range 
inhomogeneous mesoscopic structures that break the C4 
structural symmetry of the CUO2 plane. ^^"^'^ Inelastic 
neutron scattering (INS) studies of these materials re- 
veal magnetic excitations that can be described by one- 
dimensional stripe structures consistent with the STM 
results. ^^'^^ These inhomogeneous electronic structures 
may be a more general manifestations of electronic liq- 
uid crystalline phases of doped Mott insulators. ^'^^ 

Nuclear magnetic resonance (NMR) of the oxygen nu- 
clei in the cuprates offers unique insight into the inter- 
play of charge and spin degrees of freedom. ^^"^^ Unlike 
STM, NMR probes the bulk of the material, and both 
the quadrupolar moment and nuclear magnetic moment 
of the "'^^O couple to the local charge and spin environ- 
ments, respectively, ^t) (/ = |) NMR can shed light 
on the spatial correlations between the charge and spin 
degrees of freedom of these mesoscopic liquid crystalline 
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FIG. 1. (Color online) Single crystal NMR spectra of 
Bi2Sr2CaCu2 08+4 in 9 T for fields Ho applied along the 
c— axis and in the ab-plane at 110 K. Solid lines are fits as 
described in the text. 



phases. In order to investigate the local environment 
at the oxygen site, we have measured NMR spectra and 
the nuclear spin- lattice- relaxation rate (T^^) of a sin- 
gle overdoped crystal of Bi2Sr2GaGu208+5 {Tc ~ 82 
K) isotopically enriched with ^^O. This compound has 
been studied extensively via surface probes such as STM 
and angle resolved photoemission (ARPES) because as it 
cleaves easily revealing well formed surfaces normal to the 
c— direction. ^^'^■^ On the other hand, Bi2Sr2GaGu208+5 
is one of the most difficult cuprate families to investi- 
gate with NMR because a structural superlattice mod- 
ulation gives rise to a quadrupolar broadening of the 
63,65(^^j and ^°^Bi resonances. ^'^'^^ With the exception 
of a few key experiments, there have been no system- 
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atic studies of this compound with NMR. Walstedt et 
al. and later Ishida et al. reported Cu Knight shift and 
Tj~ ^ measurements of pseudogap behavior in underdoped 
and overdoped Bi2Sr2CaCu208+5 crystals;^^'^^ Takigawa 
and Mitzi reported oxygen and copper NMR evidence 
for d-wave pairing in the superconducting statc;^^ and 
recently Chen et al. reported evidence for magnetic im- 
purities based on ^^O NMR in Bi2Sr2CaCu208+5.^^ Each 
of these studies focused on the spin susceptibility mea- 
sured at the Cu and/or O sites. In this paper, in addi- 
tion to the temperature dependence of the Knight shift, 
we focus also on the temperature dependence of the elec- 
tric field gradient (EFG) and Tj~^ tensors at the planar 
oxygen site in order to investigate the possibility of C4 
symmetry breaking and the influence of the inhomoge- 
neous electronic states on the NMR properties. 



II. SPECTRAL MEASUREMENTS 

The Bi2Sr2CaCu208-i-5 single crystal was grown by a 
floating zone method. The superconducting transition 
temperature of the as-grown single crystal was 91K. Sin- 
gle crystals were enriched with ^^O by annealing for 144 
hours at 550 K in an isotopically enriched oxygen at- 
mosphere. Unlike other high Tc materials, oxygen ex- 
changes quickly in Bi2Sr2CaCu208+5 so that the ^^O 
NMR signal is detectable even in a single crystal. In 
other cuprates oxygen NMR is usually only possible in 
aligned samples, in which a powder sample (with greater 
surface area) is isotopically enriched, then mixed with 
epoxy and cured in an external field. Aligned pow- 
der samples have well resolved spectra for the applied 
field along the alignment axis (usually the c-direction), 
but exhibit a random powder pattern for perpendicular 
fields. ■^^ Although there are ways to extract parameters 
of the EFG and T^^ tensors from these powder patterns, 
there are usually unavoidable baseline and orientation 
distribution corrections. ^^'"^^ On the other hand, isotopi- 
cally enriched Bi2Sr2CaCu208+5 single crystals enable 
one to measure directly the full shift, quadrupolar and 
spin-latticc-relaxation rate tensors. This information is 
crucial because the planar oxygen does not have axial 
symmetry, and the principal axes of these tensors lie 
along the Cu-0 bond directions and the c-axis of the 
crystal. Measurements of these tensors enable one 
to extract the full anisotropy of the spin fluctuations and 
EFG at this site. 

The NMR measurements were carried out using a 
Quantum Design Physical Property Measurement Sys- 
tem (PPMS) in a field of 9 T using a home-built probe, 
and magnetization measurements were carried out us- 
ing a Quantum Design SQUID magnetometer. Spectra 
were acquired by accumulating spin echoes at fixed field, 
and the spin-lattice-rclaxation rate was measured by in- 
version recovery. The Hamiltonian for the planar 0(1) 
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FIG. 2. (Color online) The magnetic shift, (upper panel) 
and the electric field gradient (lower panel) as a function 
of temperature for a — a,b, c. The colors correspond to the 
same field orientations as in Fig. 1. The dotted line in the 
upper panel indicated Tc. 



nuclei is given by: 



u = jhHo{i+Ka)io.+-^m^-n+v{ic-m, (i) 

where 7 = 0.57719 MHz/T is the gyromagnctic ratio for 
^'Q, la and Ka are the nuclear spin operator and the 
magnetic shift along the a direction, Hq is the external 
field, i^aa are the components of the EFG tensor along 
the principal axes, and 77 = {\Vcc\ — \Vaa\)/\Vbb\ is the 
asymmetry parameter of the EFG. We adopt a conven- 
tion in which the principal axes lie along a \\ [110] in 
the Cu02 plane perpendicular to the Cu-0 bond axis, 
b II [110] is parallel to the Cu-0 bond axis (see diagram 
in Fig. 1), and c || [001] lies perpendicular to the Cu02 
plane. The largest eigenvalue of the EFG tensor (vtb) is 
associated with the bond axis. The EFG components are 
given by Vaa = SeQVaa/'^Oh, where Q = 25.6 mbarn is 
the quadrupolar moment of the ^''O, and Vaa is the sec- 
ond derivative of the potential V{r) along the a direction 
at the 0(1) site. 

Fig. 1 shows spectra of the oxygen for the field both 
parallel and perpendicular to the c axis. Note that when 
Hp lies in the plane, the two planar 0(1) sites are no 
longer equivalent, and we define 0(lb) as the site with 
the field parallel to the bond axis, and 0(la) site as 
the site with the field perpendicular. We do not ob- 
serve the oxygen site in the SrO layer because T^^ for 
these sites is much longer than the experimental rcpe- 
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FIG. 3. (Color online) Ka versus Xa for a = a,b. The colors 
match those in Fig. 1. As described in the text, the best fits to 
the data give the hyperfine coupling constants 2Caa ~ 102(40) 
kOe/^s and 2Ctt = 160(40) kOe/^s, and Koa = -0.5(4)% 
and Kob = -0.8(4)%. 



tition rate. Furthermore, structural disorder from the 
supermodulation may broaden the oxygen site in the 
BiO layer as well.^'' In order to extract the Knight shift 
and EFG tensors, we fit the spectra in each direction to 
a sum of five Gaussians. The second m oment of each 
transition is given by cr = \J <y\[ + na'q, where gm is 
the rms second moment of magnetic contribution to the 
linewidth, ctq is the rms second moment of the quadrupo- 
lar contribution to the linewidth, and n = 0,±1,±2 
indicates the particular nuclear transition. At 110 K 
we find {Ka,Kb,Kc) = (0.19(1), 0.29(1), 0.19(1))%, and 
{iyaa,nb,iycc) = (-754(3), 1139(4), -375(4)) kHz (the 
asymmetry parameter rj — 0.33). These values arc close 
to those reported previously in a crystal with similar Tc-^^ 
The EFG arises both from on-site contributions from the 
unfilled 2p orbitals of the oxygen, and lattice contribu- 
tions from distant charges. The dominant contribution 
arises from the on-site term, and is a direct measure of 
the doping in the cuprates^'''"^^. Susceptibility and resis- 
tivity measurements indicate = 82 K corresponding 
to a doping level of p « 0.21.-^^ 

Figure 2 shows the temperature dependence of Ka 
and Vaa- The EFG tensor is temperature independent, 
suggesting the absence of any static charge order. In 
contrast the magnetic shift is strongly temperature de- 
pendent. The magnetic shift is given by Ka{T) = 
K° + /iT*" + K^{T), where K°, iiT*", and are the 
contributions from the orbital (Van Vleck) susceptibility 
of the O 2p orbitals, x°; the diamagnetic susceptibility 
of the filled orbitals, xJ"j the spin susceptibility of 
the Cu 5* = 1/2 spins, XaC^)- The first two terms are 
temperature independent and combine to give an overall 
shift Kq. The spin shift is arises because of the hyperfine 



coupling to the Cu electron spins: 

^hyp=7%MB^i-C-S(r + r,), (2) 

where the sum is over the two nearest neighbor Cu sites, 
and C is the transferred hyperfine coupling tensor. "^^ The 
spin shift is given by = 2CaaX%{T), and the cou- 
pling parameters Caa can be determined by the slope 
of Ka versus Xab along the principal axes of the hyper- 
fine tensor in the normal state (see Fig. 3). We find 
Caa = 51 kOe/nB and Cbb = 80 kOe/^s, similar to 
previous measurements in YBa^CusOT^a.'^^ For the c- 
direction, Xc is essentially temperature independent, and 
we do not have enough precision to measure Ccc] how- 
ever spin-lattice-relaxation rate measurements discussed 
below suggest that Ccc ~ C'aa, in agreement with pre- 
vious reports in YBa2Cu307_5.'^^ The approximate ax- 
ial symmetry of the transferred hyperfine coupling arises 
from the O 2p-orbitals that hybridize with the nearest 
neighbor Cu 4s orbitals. '^^ 

The pseudogap is manifest in NMR as a peak in the 
Knight shift at a temperature followed by a break- 
down of scaling below a temperature T*.^^ As seen in 
Fig. 2, we find that KaiT) reaches a maximum at a 
temperature T"* ss 125 K and is suppressed at lower 
temperatures. This value agrees with other probes 
of the pseudogap temperature for this doping level in 
Bi2Sr2CaCu208+5.^''''^^ Ka continues to be suppressed 
below Tc due to the spin singlet nature of the Cooper 
pairs. 

The resonances in Fig. 1 arc broadened inhomoge- 
ncously. The fits to the spectra yield lincwidths um — 
36(8) kHz, (jQ = 98(6) kHz for the field in the plane, and 
um = 46(8) kHz, ctq = 120(6) kHz for the field along 
c. The dominant contribution to gm arises from the in- 
homogeneity of the magnet (approximately 5 x 10~^ in 
the PPMS). Fig. 4 shows the temperature dependence 
of the lincwidths of the first satellites (±| o ±\ transi- 
tion) for the field along both a and b. The linewidth is 
slightly larger for the b direction because K^ > Ka and 
therefore the inhomogeneity of the field translates into 
a greater linewidth. At low temperatures the lincwidths 
decrease slightly because K is reduced. The quadrupolar 
contribution to the linewidth, ag, reflects an inhomoge- 
neous distribution of the local EFG at the 0(1) sites. As 
discussed below, this distribution arises because of the 
intrinsic inhomogeneity in the hole occupations of the 
oxygen orbitals. 

III. RELAXATION MEASUREMENTS 

In order to investigate the spin fluctuations present at 
the 0(1) site, we have measured T^^ for all three ori- 
entations. Since 7^ for this site, the relationship be- 
tween the spin fluctuation spectrum and the spin-lattice- 
relaxation rate is complex. The eigenstates, \4>i)-, of the 
Hamiltonian (1) are superpositions of the |m) states, and 



4 



N 



CM 



0. 
0. 



to 



-8. 0. 
=o 0. 



to 



-0 



80 

60 

40 

20 


03 

02 

01 

00 

02 
00 
02 
04 



: *M 


1 
1 
1 

1 1 • H||a 
1 • H||b 

'1 1 1 1 


~ * * 1 


M 1 1 1 

" T f f I 

1 1 1 1 1 


- 


M 1 1 1 


1 


: 1 1 1 1 



50 100 150 200 250 
T(K) 



FIG. 4. (Color online) (upper panel) The average of the 
rms second moments of the first satellite transitions (| -f^ |, 
— I — ^) for the field parallel and perpendicular to the Cu- 
O bond axis, (middle panel) The rms second moment, Sup 
of the local hole concentration distribution of the two oxygen 
sites (see text for details), and (lower panel) the difference 
between Sup.a and Sup^b versus temperature. The dashed line 
indicates T* = 94 K and the dotted line indicates Tc — 82 K. 



transitions between these states are driven by fluctua- 
tions of the hyperfine field h.{t). This field arises through 
the interaction in Eq. 2, and gives rise to a fiuctuating 
Hamiltonian: 



Hiit) ^ jh ha{t)Ia + hb{t)h + hc{t)Ic 



(3) 



where h(t) ~ gfJ-sC ■ J^i^i^ + ^i,t)- The spin fiuc- 
tuations drive the nuclear spin lattice relaxation rates 
Ra = {ha{t)ha{t + T))e~'^'^'^dT and we assume 

that h(i) fiuctuates independently in each direction. The 
detailed form of the magnetization relaxation and the de- 
pendence of T^^ on Ra are derived in the appendix. The 
measured values are shown in Fig. 5 for all three direc- 
tions. The relaxation rates in each direction scale with 
one another and are reduced at low temperature due to 
the pseudogap. As seen in Fig. 5(b), the fiuctuations 
are largest along b, in agreement with the anisotropy of 
the hyperfine coupling. Fig. 5(c) displays Ra/T and 
{K^)^ for each direction. The scaling between these two 
quantities reflects Korringa relaxation {TiT ~ (/ST*)^), in 
which the nuclei are relaxed by spin-flip scattering with 
quasiparticles, and agrees with previous studies. ^^'^^'^^'■^^ 
Fig. 7 displays (Ra/T)^^'^ versus K with temperature as 
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5. (Color online) (a) the spin-lattice-relaxation rate, 
T^^, for the field along each direction a (•, red), b (■, blue), 
and c (a, green), (b) The relaxation rates Ra determined 
from the Tia data as described in the text, (c) Ra/T and 
{K^f (solid lines) versus T. Here = K - Kq, where Kq 
is determined in Fig. 7. 



an implicit parameter. The fact that the data are linear 
further supports this interpretation and implies that the 
temperature dependence of T^^ of the 0(1) site is driven 
by a partial gapping of the density of states (pseudogap) 
rather than by spin fluctuations. 



IV. DISCUSSION 

A. Scaling behavior 

In recent years evidence has emerged that the the mag- 
netic behavior of the YBa2Cu307_5 and La2-£i;Sra;Cu04 
cuprate families is best understood in terms of two elec- 
tronic degrees of freedom. Over a broad range of 
doping, the uniform susceptibility scales as x(^7 ^) = 
Xo{x) + x"x(T/T™(x)), where x(T/r"(x)) is a univer- 
sal function of T/Tm, and T™ is the temperature where 
X is maximum. In Fig. 6 we show the Knight shift scaled 
to fit this form for all three directions using T™ = 125 
K. The data scale well for T/T™ > 0.75, suggesting 
that T* ^ 94 K. Barzykin and Pines have found that 
T" = 1218ii:(l - 4.45p), where p is the doping level in 
La2-_a:Sr3.Cu04. Using the measured value of T"* we find 
p = 0.201, which is close to the value of 0.21 based on 
the measured Tc- These results confirm that the mag- 
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netic scaling in Bi2Sr2CaCu208-)-5 agrees quantitatively 
with that of La2_a;Sr2.Cu04, and implies that the mag- 
netic scaling is uniform across all high Tc families. 



B. Evidence for electronic inhomogeneity 

In recent years an increasing number of experiments 
on Bi2Sr2CaCu208+5 and other high temperature su- 
perconductors have uncovered evidence that the elec- 
tronic degrees of freedom develop mcsoscopic ID stripe- 
like structures of inhomogeneous spin and charge densi- 
ties below the pseudogap energy/'^^'^'^^^^'^^''*^ STM mea- 
surements have revealed that the 90° C4 rotation sym- 
metry of the Cu02 plane is broken by the electronic 
degrees of freedom in underdoped Bi2Sr2CaCu208+5 at 
the pseudogap energy, and that the two crystallographi- 
cally equivalent planar O sites acquire different local hole 
concentrations.*'^ In this case the local EFG and Knight 
shift should differ for the two sites, leading to a split- 
ting or broadening of the NMR resonances. These asym- 
metric patterns are short-ranged, may be dynamic, and 
may or may not couple to the external magnetic field 
used for the NMR experiments. In order to investigate 
any asymmetry in the NMR response of the two planar 
O sites, we consider the linewidths of the parallel and 
perpendicular sites. There are three potential contribu- 
tions to the linewidth: (i) inhomogeneous fields from the 
magnet, (m) a distribution of local Knight shifts, and 
(iii) a distribution of EFGs due to locally varying hole 
concentrations in the O 2p orbitals. The field inhomo- 
geneity /3h = SH/H « 5 X 10"'* is independent of di- 
rection and temperature, whereas the contribution from 
the Knight shift /3_r- = SK/K depends on both quanti- 
ties. The distribution of local EFGs is directly related to 
the hole occupation in the O 2p orbitals via the relation 
Uaa = Vaa + <Za'^■p^ whcrc ly^^ are material dependent 
constants, qa = 1.225 MHz for a = (a,c), and 2.452 
MHz for a = b, and Hp is the hole concentration in the 
O 2p orbital.'^* Assuming the holes go exclusively into 
each of these two orbitals, Up = p/2 = 0.105. For a dis- 
tribution Sup of local hole doping, the quadrupolar con- 
tribution to the linewidth would be ctq ~ nqaSup, whcrc 
n is the particular satellite transition. If each of these 
contributions to the line broadening is independent of 
one another (which may not necessarily be the case (see 
Ref. 18)), then we can estimate the rms second moment 
of hole concentrations in both the parallel {a = b) and 
perpendicular (a = a) O 2p orbitals as: 



Sn 



■p,a{T) — 

nqa 



^al{T)~ul{Pl+PlKl{T)), (4) 



where vq = 7i?o- By fitting the central transition (for 
which (Tq = 0) at llOK we estimate Pk = 0.0024 for 
a = a and 0.0036 for a = b. We then calculate the tem- 
perature dependence of Snp{T) for both directions, as 
shown in the middle panel of Fig. 4. This result shows 
that the distribution of local hole dopings is roughly equal 
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FIG. 6. (Color online) The spin component of the Knight 
shift versus T/T", where = 125 K. The solid line is 
the Johnston-Nakano scaling form'^'*. 



to 10-15% of the average value, in agreement with pre- 
vious STM estimates, and confirming that this inhomo- 
geneity observed by STM is also present in the bulk.^^'^^ 

On the other hand, we find no evidence for any nematic 
asymmetry between the hole distribution for the parallel 
and perpendicular sites (lower panel in Fig. 4). Note, 
however, that the STM nematic asymmetry was observed 
in zero external field at finite bias in heavily underdoped 
samples. This NMR study was conducted in a field of 
9 T in a slightly overdoped sample, and any dynamic 
nematic fiuctuations would likely give rise to motional 
narrowing of the resonance. Furthermore, it is unclear 
what role the magnetic field would have on the pattern 
of electronic excitations. Further studies in underdoped 
samples may help clarify these issues. 

Surprisingly, there appears to be little or no contribu- 
tion to the linewidth from the structural modulation.^* 
This modulation dramatically broadens the Cu and 
Bi resonances, as well as the oxygen site in the BiO 
layer. ^^'^^'** The Cu02 planes are buffered from the dis- 
ordered BiO layer by the SrO layers, and the relatively 
small quadrupolar moment of the oxygen render this site 
less sensitive to the structural modulation. 



C. Korringa relaxation 

The scaling between the Knight shift and the spin lat- 
tice relaxation rates Ra/T observed in Figs. 5 and 7 
is striking because it implies that the spin fluctuations 
responsible for the oxygen relaxation arise from quasi- 
particle scattering on the Fermi surface, rather than from 
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FIG. 7. (Color online) ^RajT versus Ka for a = a, 6, c. The 
colors match those in Fig. 1. The solid line is the best fit to 
the data as described in the text with slope 103(11) (s K)^''^. 
The dotted line has a slope given by the Korringa value (see 
text). 



fluctuations of the Cu spins. In general, the nuclear spin- 
lattice-relaxation rate is given by: 

where x" ('^j q) is the imaginary part of the electron spin 
susceptibility.^^'^^ The hyperfine form factor FQ,(q) = 
2Caa cos((72-,ya/2) filters out Cu spin fluctuations at 
Q = (tt/o, 7r/a) because the oxygen is located symmetri- 
cally between the two coppers. ^^'"^^ However recent neu- 
tron scattering studies on Bi2Sr2CaCu208+5 have con- 
firmed the presence of incommensurate spin excitations 
at wavevectors Qo 7^ Q.^^ In this case the ^^(Qo) 
does not vanish, and spin fluctuations should contribute 
to the spin lattice relaxation rate. The temperature 
dependence of i?Q would be driven by the correlation 
length, ^(T), and Ra/T would not scale with the shift 
Ka ^ x{'^ = 0, q = 0).'^^''*'' This scenario is inconsistent 
with the data. 

On the other hand, for relaxation driven by quasiparti- 
cle scattering at the Fermi surface, ^/RajT ~ Ka/V^, 
where k = ^^/Trksfvy'^ is the Korringa constant.^* The 
best fit in Fig. 7 yields a slope of 103(11) (s K)-'^/'^, which 
is remarkably close to the theoretical value 187 (s K)"-'^/^, 
suggesting that quasiparticle scattering rather than an- 
tiferromagnetic fluctuations is the dominant mechanism 
for the oxygen spin lattice relaxation. This result sug- 
gests that either (1) the oxygen is completely insensi- 
tive to the spin fluctuations, (2) there is a second elec- 
tronic degree of freedom that couples to the oxygen, or 
(3) there or no antiferromagnetic spin fluctuations in the 



overdoped regime. Since ^ '--^ a at T = T,„, it is likely 
that short range spin fluctuations continue to exist down 
to TJ' 

The discrepancy between the neutron scattering and 
the oxygen NMR data is well known, and others have 
proposed several possible explanations, including (1) the 
presence of another band or degree of freedom associated 
with the holes on the oxygen orbitals,^^''^'^'^^ (2) a sec- 
ond transferred hyperfine coupling to the next nearest Cu 
sites,"'^ and (3) fluctuations of the EFG at the O site.'^" 
Another possibility is that the dynamical susceptibility 
contains two terms, one from the Fermi surface itself, 
and another from the incommensurate antiferromagnetic 
background. "^^ The incommensuration response implies a 
modulation in real space of the staggered magnetization 
with wavelength A ^ a.^^ If the response were static, 
then there would be multiple oxygen sites with different 
hyperfine fields and relaxation rates. The experiments 
clearly rule this static scenario out, thus we expect dy- 
namic fluctuations of this modulated structure. In this 
case the oxygen sites would be motionally narrowed, and 
all of the oxygen sites would experience antiferromagnetic 
fluctuations that would drive the spin lattice relaxation. 

If this modulation contained higher order harmon- 
ics such that the real space spin density were a square 
wave rather than a sinusoid (discommensurations) the 
response would be commensurate with antiphase domain 
walls. ■^^'"'^ In this case the oxygen form factor would filter 
out the spin fluctuations, which resolves the apparent in- 
consistency between the INS and NMR results. In this 
case the relaxation of the oxygen may be driven by the 
dynamics of the domain walls. An alternative explana- 
tion is that the relaxation is driven by a second electronic 
degree of freedom. Recent NMR studies of the O and Cu 
shift in La2-a:Sr3:Cu04 suggest the presence of a second 
spin component at the oxygen site."^^'^" Which explana- 
tion hold for this case remains unclear at present. 



V. SUMMARY 

We have measured the full anisotropy of the EFG, 
Knight shift, and spin lattice relaxation tensors at the 
planar oxygen site in an overdoped single crystal of 
Bi2Sr2CaCu208+5. We present a complete analysis of 
the anisotropy of the spin fluctuations on the spin lattice 
relaxation rate for the oxygen. The Knight shift and spin- 
lattice-relaxation rate reveal the presence of a pseudogap 
and the relaxation rate is well described by quasiparticle 
scattering at the Fermi surface. The Knight shift agrees 
with the scaling result implying that Tm ^ 125 K and 
T* ~ 94 K. Antiferromagnetic fluctuations of Cu spins do 
not appear to contribute to the oxygen relaxation rate, 
suggesting that the incommensurate order observed by 
neutron scattering must contain higher order harmonics. 
The quadrupolar contribution to the resonance linewidth 
is consistent with a distribution of local hole concentra- 
tions in the oxygen 2p orbitals, in agreement with STM 
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studies. This result implies that the intrinsic electronic 
inhomogeneity is present in the bulk and is not purely 
a surface phenomenon. On the other hand, the charge 
distribution in the O 2p orbitals appears to be identical 
for both the planar oxygen sites and remains tempera- 
ture independent. Further studies at lower dopings in the 
pseudogap phase may provide new information about the 
interplay of the local spin structure and the complex elec- 
tronic patterns that emerge in the underdoped cupratcs. 

Appendix 

Spin lattice relaxation is described by the master equa- 
tion: 

^ = iy,,(p,,(t)-p„0), (A.l) 

where pu are the diagonal components of the thermal av- 
eraged density matrix, is the thermal equilibrium den- 
sity matrix, and Wij is the spin lattice relaxation matrix. 
The elements of Wij are given by: 

where 

/oo 
{c^^\nl{t)\(j)J){c^l\nl{t + r)|(/)fc)e-'""dr 
-OO 

(A.3) 

and the overline signifies a thermal average.'*® If [IzjH] = 
such that l^,) are eigenstates of I^, then the Jijki van- 
ish for any transition where j ^ i ± 1 and fc ^ ^ ± 1 as 



per the usual selection rule. For the 0(1) site the Jijki 
do not vanish because the eigenstates are superpositions 
of the states. Effectively this means that for the field 
along a direction a, fluctuations of h{t) along all three di- 
rections can give rise to transitions instead of just those 
two directions that are perpendicular to a. The relax- 
ation matrices, Wa, depend simultaneously on Ra, Rb, 
and Re- This simultaneous dependence means that the 
coupled differential equations in (A.l) cannot be writ- 
ten in terms of a single effective time constant (Ti), and 
there is no closed form expression for the magnetization 
recovery function. 

In order to surmount this difficulty, we expand the re- 
laxation matrix as: 

Wa{Ra, Rb, Rc) — T^^aiRa, Rb, Rc)Wa{Ra/Rb, R^al Rc) + 

SW^iRa,Rb,Rc), (A.4) 

where WaiR^/Rb, Rc/Rb) ^ constant matrix, 

T^^{Ra, Rb, Rc) is an overall time constant that depends 

on the relaxation components Ra, and SWa{Ra, Rb,Rc) 
is a small correction which we ignore. We estimate the 
ratios = 0.329 and = 0.373 as the squares 

of the hyperfine couphng ratios, since i?Q ~ /i^ ~ C'aa: 
where Caa is the hyperfine coupling. '^^'^^ In this approx- 
imation, we assume a single component coupling model 
to isotropic Cu S" = 1/2 spins. In the limit jHo/i/bb oo 
the eigenstates of (1) indeed commute with 1^ and the 
relaxation matrices reduce to their conventional forms 
that depend only on a single time constant. In our case 
jHo/h'bb ~ 50 so the correction term is small and the 
approximation is valid. In this case we find 



/ -2.452 

-2.497 

2.497 

2.452 

0.0001201 

V 0.0001099 



2.452 

2.497 

-6.472 

-6.392 

3.975 0.0002036 

0.0002096 3.940 



0.0001099 \ 

0.0001201 

3.975 0.0002096 

0.0002036 3.940 

-8.427 4.451 

4.451 -8.391 / 



Wb 



and 



/ -2 



V 0.00004005 





-2.487 
2.487 


0.00004574 




Wc 



-2.440 



2.440 


V 0.0002059 



2.488 0.00004005 \ 

2.487 0.00004574 

-6.466 3.979 0.00007869 

-6.469 0.00007528 3.980 

3.979 0.00007528 -8.456 4.477 



0.00007869 3.980 





-2.501 
2.501 


0.0002153 






2.501 
-6.477 


3.975 
0.0003816 



2.440 



-6.366 
0.0003760 
3.926 



4.477 





0.0002153 

3.975 
0.0003760 
-8.420 
4.444 



0.0002059 \ 


0.0003816 
3.926 
4.444 
-8.371 / 



-8.458 / 



(A.6) 



(A.7) 
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where 
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The magnetization recovery measured along the a 
direction following an inversion pulse is then given 
by: Ma{t) = Mo(l - 2fMt/Tia,)), where Mx) = 
X]i=i : and the constants Ai and /3i are given in 

Table I for the first satellite transition of the 0(1). 



TABLE I. Coefficients and exponents in the multiexponential 
spin lattice relaxation decay function for the 0(1) nucleus. 



i 


A. 


13a 


Xt 


/36 






1 


0.454 


14.84 


0.446 


14.92 


0.457 


14.82 


2 


0.440 


9.89 


0.447 


9.95 


0.438 


9.88 


3 


0.0225 


5.94 


0.0251 


5.97 


0.0215 


5.93 


4 


0.0546 


2.97 


0.0535 


2.98 


0.0549 


2.96 


5 


0.0282 


0.99 


0.0286 


0.99 


0.0280 


0.99 
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